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ABSTRACT 

We present an extensive set of stellar population synthesis models, which self-consistently 
include the (optical-far-UV) continuum emission from stars as well as the resulting emission-line 
spectrum from photoionized gas surrounding massive stars during their main- sequence-life time. 
The models are presented as a compiled library, ADEMIS, available electronically. ADEMIS contains 
the equivalent widths and the intensities of the hues [O II] AA3727, H/3, [O III] A5007, Ha, and 
[N II] A6584, which were calculated assuming a metallicity of 0.2, 0.4, 1.0, and 1.5 Zq and a wide 
range of ionization parameters. We investigate the regime of continuous star formation, assuming 
a Salpeter initial mass function, whose upper mass limit is an input parameter. The calculated 
equivalent width models, which are function of the burst age, are compared with the Jansen et al. 
atlas of integrated spectra for nearby galaxies. We reproduce the observed properties of galaxies 
along the full Hubble sequence and suggest how the metallicity and age of such galaxies might 
be roughly estimated. 

Subject headings: galaxies:general - galaxies: star forming - galaxies: emission line - ISM: HII regions 

1. Introduction newly formed stars versus the 'excited' interstellar 

medium (ISM), in a single observable quantity, the 

The presence of emission lines is a conspicu- integrated optical spectrum. It is for this reason 

ous feature of optical spectra of any galaxy that ^^iiit the modeling of the emission-line spectra oc- 

is currently forming stars. The lines result from curring in star-forming galaxies must include the 

the reprocessing of the UV ionizing radiation from interaction between both phases (stars + gas) in 

the most massive of the young stars. Therefore, ^ self-consistent manner. Partly for this reason, 

emission lines carry imprints of both phases, the ^his paper will emphasize the study of equivalent 



widths"'^ (EWs or Wx). The optical stellar contin- 
uum underlying the emission lines consists of much 
more than just the stars that ionize the gas. Actu- 
ally, in any evolved stellar population, the optical 
and red continua are dominated by red giants and 
low-mass main sequence stars. In order to model 
the observed continuum properly, one must there- 
fore consider the past stellar formation history of 
the whole stellar population, an essential step if 
the emission-line equivalent widths are to be cal- 
culated reliably. 

Working toward these goals, Garcia Vargas, 
Bressan & Di'az (1995), Stasinska & Lcithcrer 
1996, Stasinska, Schaerer & Leitherer (2001) built 
models in which the spectral energy distributions 
(SEDs) generated by stellar evolutionary synthesis 
codes were used as input to photoionization calcu- 
lations of emission-line strengths. These models 
aimed mainly at constraining the ages and mctal- 
licities of isolated extragalactic H II regions that 
are characterized by the absence or weakness of 
an old stellar population, suggesting a prosont-day 
strongly enhanced star formation rate. In short, 
these earlier models followed the evolution of spec- 
tral properties of an instantaneous burst of star 
formation. Recently, Kewley et al. (2001), Zack- 
risson et al. (2001) and Moy, Rocca-Volmerange & 
Fioc (2001) introduced models in which the star 
formation rate is constant for a few megayears, 
without interruption. These were used to derive 
the properties of star burst galaxies. 

In an attempt to include the forgotten dust 
phase, Magris (1993) combined the spectral evo- 
lution synthesis code of Bruzual & Chariot (1993) 
with the photoionization code MAPPINGS i (Bi- 
nette, Dopita & Tuohy 1985). The particularity 
of these models is that they included the effects of 
dust on the line- and continuum-omission in a self- 
consistent manner. On the other hand, only the 
solar metallicity case was studied, since detailed 
non-solar mcrtallicity tracks and stellar spectral li- 
braries wore not complete at that time. 



Equivalent widths present the advantage of relating any 
line flux with an observable quantity at the same wave- 
length, namely the underlying continuum. There are quite 
a few papers in the literature that are concerned with the 
calibration of the [O II] equivalent width and whether it 
can be used as a reliable tracer of star formation. Such 
diagnostics of W;i([0 II]) in the UV would be particularly 
valuable at redshifts where li/3 or Ha become inaccessible. 



Barbaro & Poggianti (1997) have presented a 
model that is used to calculate the emission lines 
of [O II] AA 3727 and of the Balmer series along 
with photometric quantities. The purpose was to 
estimate the present day star formation rate in 
late type galaxies. In their model, the star forma- 
tion rate and the chemical evolution of the stellar 
phase for each Hubble type have been derived from 
their own chemical evolution model, assuming, on 
the other hand, a fixed solar composition for the 
ionized gas. 

Chariot & Longhetti (2001) computed the line 
and continuum emission from star-forming galax- 
ies and presented a family of estimators of the star 
formation rate, the nebular oxygen abundance, 
and the optical depth of the dust, in galaxies, as 
a function of the available spectral information. 
Their model is generally successful in reproducing 
the observed properties^ of template galaxies of 
different types. Their analysis, however, did not 
include the modeling of equivalent widths that are 
particularly suitable for analyzing the interrelated 
properties of the stellar and nebular phases. 

We have developed a model along similar lines, 
which calculates, for an arbitrary star formation 
rate, both the nebular emission-line spectrum of 
the H II regions associated to the newly born mas- 
sive stars and the SED of the underlying stellar 
population. However, rather than discussing line 
ratios and stellar energy distributions separately, 
our work will emphasize the relationship between 
the two phenomena, by presenting our model re- 
sults as a library, ADEMIS, of emission-line equiv- 
alent widths. As supplementary information, the 
tables also include the line luminosities of the most 
prominent species observed in normal star- forming 
nearby galaxies, i.e. [O II] AA3727, H/3, [O III] 
A5007, Hq, and [N II] A6584. The models are char- 
acterized by a metallicity range of 0.2, 0.4, 1, and 
1.5 times solar (gas -|- star) and consider rates of 
star formation that are constrained by both the 
line and the stellar continuum emission, of galax- 
ies of all Hubble types. 

The method used to derive our evolutionary 
models is presented in §2. A description of the 
results from ADEMIS follows in §3 accompanied by 
a comparison of the models with the Jansen et al. 

■^Scc Chariot ct al. (2002) for a more recent application of 
this model to the determination of the star formation rate. 
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(2000a,b; hereafter JFFC) atlas of nearby galax- 
ies. The conclusions are summarized in §4. A 
more extensive version than the adjoining tables 
of ADEMis models is available on the Web. 

2. Computational Method 

The method used to derive the nebular syn- 
thetic spectrum consists first in obtaining the ion- 
izing spectrum corresponding to the hot main - 
sequence stars of an arbitrary stellar population. 
This ionizing spectrum is subsequently used by a 
photoionization code to compute the emission-line 
spectrum. The line fluxes are then combined with 
the integrated continuum spectrum of the stel- 
lar component to provide us with the equivalent 
widths of the emission lines of interest. 

2.1. The Population Synthesis Code 

One of our aims is to derive a nebular spec- 
trum fully consistent with the young stellar pop- 
ulation from which arise H II regions within an 
evolved normal galaxy. Special care must there- 
fore be given to the population synthesis of the 
star-forming regions. For this purpose, we use the 
latest version of the evolutionary population syn- 
thesis code of Bruzual & Chariot (2003, hereafter 
BC03). 

With this code, we can derive the spectropho- 
tometric properties of any specified stellar popu- 
lation. For a given star formation rate (SFR) and 
initial mass function (IMF), a model consists of a 
set of spectra of the integrated stellar population 
calculated at 221 different ages t. Each such spec- 
trum is defined at a common set of wavelength 
points covering the range 5 100 /im. The BC03 
code, furthermore, can store separately the contri- 
bution of stars in different evolutionary stages. In 
the current work, it was necessary to generate an 
integrated spectrum for all the stars and another 
one for the subpopulation of main-sequence (MS) 
stars only. This latter spectrum was essential to 
derive the ionizing SED, which excites the H II re- 
gions and whose contributors consist exclusively of 
massive MS stars. This procedure excludes, for in- 
stance, other types of hot stars that generate ion- 
izing photons, such as planetary nebulae nuclei, 
white dwarfs, or core helium-burning stars, which 
are irrelevant to young H II regions. In a separate 
paper, we have verified that the small Ha widths 



of normal ellipticals might be explained by pho- 
toionization from post-AGB stars alone (Binette 
et al. 1994). 

The time-evolution of the stellar population is 
followed using the Girardi et al. (2000) stellar 
tracks for low- and intermediate-mass stars (i.e. 
from 0.15 to 7 Mq), which are the most recent 
set from the Padova group, and the Bertelli et al. 
(1994) tracks for stars from 9 to 120 Mq. The 
advantage of these models with respect to ear- 
lier ones is that the data-base is computed with a 
homogenous input physics for all stellar tracks in 
the grid. For instance, moderate convective over- 
shooting from convective cores and mass loss from 
massive stars are taken into account. In the most 
recent version used here, these tracks have been 
improved further with the adoption of new low- 
temperature opacities (see Girardi et al. 2000 for 
details) 

The Kurucz stellar atmospheres spectra library 
(distributed by Lejeune, Cuisinier, & Buser 1997) 
has been adopted to compute the stellar proper- 
ties of stars with Tefj < 50, 000 K. For the higher 
temperature stars, in order to give an adequate 
treatment to the ionizing continuum of hot MS 
stars, we use the Ranch (1977) NLTE model at- 
mospheres'^ for A < 2000 A, and a simple black 
body spectrum for larger wavelengths. We rec- 
ognize that large differences exist in the ionizing 
flux between LTE model atmospheres that include 
only H, He, C, N, O, and Ne, (like the Hummer 
& Mihalas (1970) models used in the emission-line 
synthesis of Magris (1993)) and the NLTE Rauch 
(1997) models, which include the opacity of heav- 
ier elements from Li to Ca. Despite these differ- 
ences, we find, however, that the general prop- 
erties of the ionizing flux that emerges from MS 
stars in a composite stellar population are not sig- 
nificantly affected. The reason is that high-mass 
stars have a relatively short lifetime for Tg// hot- 
ter than 50, 000 K, and their contribution to the 
integrated spectrum of a composite stellar pop- 
ulation becomes insignificant for any lasting star 
formation burst. 

2.2. The Photoionization Code 

To compute the nebular spectrum, we use 
the multipurpose shock/photoionization code 

^Also available at http://astro.uni-tuebingen.de/~rauch/. 
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MAPPINGS ic described by Fcrruit ct al. (1997); 
see also Binette et al. 1985), in which the ion- 
ization and thermal structure of the nebula are 
computed in a standard fashion, assuming local 
ionization and thermal equilibria throughout the 
nebula. The calculation is made at a progression 
of depths in the nebula, until all the ionization 
radiation has been absorbed. 

For the nebular gas distribution, we adopt a 
simple slab geometry as in the work of Shields 
(1974), who argued that the integrated line spec- 
trum was not significantly different from that cal- 
culated with the more usual spherical geometry. It 
is also the geometry preferred for the Orion neb- 
ula as modeled by Baldwin et al. (1991) or Bi- 
nette, Gonzlez-Gmez & Mayya (2002). The ratio- 
nale is that H II regions are rather chaotic in ap- 
pearance, showing evidence of strong density fluc- 
tuations with much of the emission coming from 
condensations or ionized "skins" of condensations. 
The relation between the dimensionless ionization 
parameter U, as defined in this work*, and the 
emission lines, is also more straightforward and 
intuitive in the plane-parallel slab case. 

Our set of solar abundances is taken from An- 
ders & Grevesse (1989). When varying metallic- 
ity, we will assume that the abundances scale with 
metallicity in a linear fashion for all elements ex- 
cept helium and nitrogen. The He/H ratio will 
be kept constant at the solar value. The observed 
behavior of the N/0 ratio in local and external 
H II regions suggests that the nitrogen abundance 
increases via secondary nucleosynthesis for metal- 
licities above 0.23 solar (van Zee, Salzer & Haynes 
1998; see Henry and Worthey 1999 for a review). 
To represent such N/0 behavior, we adopt the fol- 
lowing relations (Dopita et al. 2000): 

log(N/H) = -4.57+log(Zg„,/ZQ); \og{Zgas/ZQ) < 

log(N/H) = -3.94 + 2 \og{Zgas/ZQ)- otherwise 
where Zgas denotes the gas phase metallicity. 

^Wc define the ionization parameter as follows, U = 
Q£fo/(47rr^cn£f ), where Q^o is the photon luminosity 
(quanta s~^) of the stellar cluster, c is the speed of light, 
and r is the distance separating the photoionized gas of 
density nn and the ionizing stellar cluster. 



2.3. The Jansen et al. atlas as comparison 

sample 

We will compare the predictions of our model 
with the spectrophotometric atlas of nearby galax- 
ies of JFFC, which is ideal for this purpose, since it 
includes a representative sample of galaxies span- 
ning the whole Hubble sequence and extending 
to lower luminosity systems than the Kennicutt 
(1992) sample. This sample is frequently referred 
to as the Nearby Field Galaxy Survey (NFGS). 
Such a comparison with the ,1FFC data will al- 
low us to test the global properties of our models 
in the observational framework of a wide range of 
galaxy luminosities. The full sample includes 196 
nearby galaxies from E to Irr, of luminosities in 
the range Mb = -14 to -22. However, of the 141 
objects with detected Ha emission, we will con- 
sider only those galaxies for which the equivalent 
widths EW(Ha) > 10 A and EW(H/3) > 5 A. 
This will ensure a reliable measurement of Balmer 
decrements according to Jansen, Franx, & Fab- 
ricant (2001) and leave us with a sample of 85 
objects. The atlas contains the equivalent widths 
and line ratios with respect to H/3 of the following 
lines: [O II] AA 3727, RS, H/3, [O III] A4959, [O III] 
A5007, He I A5876, [O I] A6300, [N II] A6548, Ha, 
[N II] A6584, [S II] A6716, and [S II] A6731. 

Since Jansen et al. (2001) have reported that 
the Balmer emission-line fluxes and equivalent 
widths given in JFFC were only partially corrected 
for stellar absorption, we have applied the addi- 
tional correction suggested by these authors of 1 
and 1.5 A (EW) for H/3 and Ha, respectively. 

We corrected all line ratios for interstellar red- 
dening by comparing the observed Balmer decre- 
ment, Ha/H/3, with the case B recombination 
value of 2.85 (appropriate for Ue ^ 10^ — 10'*cm~^ 
and T = lO^K: Osterbrock 1989), and using the 
,^andy et al. (1975) Galactic extinction curve as 
tabulated by Seaton (1979). 

2.4. Grid of photoionization calculations 

Photoionization models were calculated using 
the various ionizing SEDs computed by the BC03 
code. These SEDs differ by the age t of their 
stellar population and correspond to MS stars, 
which were all formed at t = 0, following a sim- 
ple Salpeter (1955) IMF that covers a mass range 
comprised between a lower and upper mass limit, 
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miow and m„p, respectively. The adopted lower 
limit is miow = 0.09 M© in all models. Because of 
the crucial role played by the upper mass limit on 
the hardness of the energy distribution, we have 
explored a wide range of upper mass limits: m„j, = 
120, 100, 90, 80, 70, 60, 50, and 40 M©. As for the 
mctallicity of the stellar population, , we have 
calculated solar metallicity models (Z* = 0.019) 
as well as metal-poor (Z* = 0.008 and 0.004) and 
metal-rich (Z* = 0.03) models. 

2.4-1- Usefullness of the hardness parameter t] 

We measure the hardness of the ionizing spec- 
trum using the ratio between the flux of stellar 
photons (quanta s^^) emitted above 24.6 eV and 
those emitted above 13.6 eV, jj = QHe°/QH°- 
The hardness rj can be associated with the cluster 
equivalent effective temperature^ frequently used 
by H II region modelers (e.g. Shields & Tinsley 
1976; Stasinska 1980; Campbell 1988). 

In Fig. la, we show the hardness, ry, as a func- 
tion of rUup, of various ionizing distributions that 
arc later used in photoionization calculations. An 
instantaneous burst with a Salpeter IMF was as- 
sumed. Each curve corresponds to a different stel- 
lar track metallicity, Z*, as labeled. The differ- 
ent curves show a reduction in rf with increas- 
ing metallicity, because the higher Z^,, the cooler 
the stellar evolutionary tracks, and therefore, the 
softer the SED (and the lower 77). This behavior 
has previously been reported by McGaugh (1991), 
who used the zero- age point of the stellar evolution 
tracks of Maeder (1990) to compute the ionizing 
radiation field of a cluster of stars distributed as 
the Salpeter (1955) IMF. Even though the details 
about our respective ionizing spectra may differ, 
it is found in this work as in McGaugh, that for 
a given SED with a certain 77, it is possible to 
encounter another SED with the same 77 if we de- 
crease (or increase) appropriately both rUup and 
Z,:. This interesting property, inferred from in- 
stantaneous burst models, tells us that both pa- 
rameters are in fact degenerate. 

Figure lb shows the temporal evolution of 77 for 
a model with a continuous SFR calculated with 
a Salpeter (1955) IMF, m„p = 125M0, and the 
metallicities as labeled. We can see that rj reaches 
a constant asymptotic value after ss 10^ yr, due to 
the equilibrium between the number of new born 
massive stars and the number of them evolving off 



the MS. The same kind of evolution occurs with 
any star formation rate provided it is continuous 
and either constant with time or monotonically 
decreasing. In Fig. lb, we also show the behavior 
of r]{t) for the case of a solar metallicity instanta- 
neous burst, for comparison (dashed line). 

In our extensive grid of photoionization cal- 
culations, U was varied over the wide range of 
—4 < log U < —1.5 although, in practice, a more 
limited range suffices, as independently shown in 
the modeling of H II regions by different authors 
(e.g., Evans & Dopita 1985; Campbell 1988). 

As for the nebular gas metallicity, Zgas, as a 
rule we adopt the same value as for the stellar 
atmospheres, and we will thereafter use Z to refer 
to the metallicity of either component. 

As a result of these computations, wc obtained 
an extensive grid of emission-line spectra parame- 
terized according to U, r], and Z. We verified that 
it was sufficient to read the appropriate emission 
line fiuxes from our grid of photoionization mod- 
els, using T] as the lookup value. Models at in- 
termediate values of the parameters U , rj, and Z 
are obtained by simple interpolation (in the log- 
plane) within the grid points. This procedure al- 
lowed us to do away with calculating separate new 
photoionization calculations each time r] changed. 

2.4-2. Two standard diagnostic diagrams 

In Figure 2 we show a standard diagnostic dia- 
gram of the relationship between 
[O IIl]A5007/H/3 and [O II]AA3727/H/?. The solid 
lines join discrete models (open circles) calculated 
by using the rUup value that appears in panel b 
near the right extremity of each line. Different 
values correspond to different 77- values with 77 
strongly dependent as well on metallicity as Fig- 
ure 1 already showed. The dotted lines join points 
of constant U (indicated as log ?7-values near the 
upper end of the curves). With a range of ioniza- 
tion parameter of about a factor of 10 (—3.25 < log 
U < —2.25) as plotted in this figure, and a range of 
metallicity between 0.4 and 1.5 Zq, our photoion- 
ization calculations encompass the dispersion in 
the [O III]/H/? versus [O II]/H/3 relation observed 
in the JFFC sample of galaxies. 

The same range in ionization parameter and 
mctallicity characterizes the cxtragalactic H II re- 
gion sequence of Dopita et al. (2000). Similarly, 
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Bresolin, Kcnnicutt & Garnctt (1999), who used 
a diflFerent but scalable definition of U, applicable 
to their geometry, concluded that a range of only 
a factor of 10 in ?7 suffices to reproduce the opti- 
cal emission-line spectrum of their sample of spiral 
galaxy H II regions. Our models also happen to 
overlay the correlation presented by Kobulnicky, 
Kennicutt, & Pizagno (1999) in their [O III]/H/? 
versus [O II]/H/3 diagram, which superposes indif- 
ferently individual H II regions and galaxy inte- 
grated ratios. 

In Figure 3 we show the line ratio diagnostic 
diagram of [O III]/[0 II] versus [N II] /[O II], which 
Dopita et al. (2000) suggests as being optimal for 
separating variations in the ionization parameter 
from variations in chemical abundance. The upper 
panel has been calculated with a simple zero-age 
stellar population, while the lower panel used a lO'' 
y old population formed according to a continu- 
ous SFR prescription. The [O III]/[0 II] ratio has 
been established to be useful for determining the 
ionization parameter (e.g. McGaugh 1991), while 
the [N II] /[O II] ratio was found to be virtually in- 
dependent of the ionization parameter, albeit very 
sensitive to the abundance for metallicities above 
0.23 solar. This Figure indicates that the galaxies 
in the JFFC sample span a range in abundances 
from 0.2 solar up to solar or possibly 1.5 times so- 
lar, without any evidence for much higher metal- 
licity objects^. As in Figure 2, log U spans the 
range —3.25 to —2.25. 

Interestingly, the continuous SRF models of 
Fig. 3b overlap better the observational points 
than the zero-age models of Fig. 3a. This simply 
results from having a softer SED (smaller rf) in the 
continuous SFR case, which produces models with 
a lower [O IIl]/[0 II] ratio. 

2.4-3. The synthesis of line spectra, stellar spec- 
tra and EWs 

The SED of a given stellar population with 
age t and arbitrary SFR, IMF, and Z is com- 

^Thc alternative [N II]/Ha ratio has previously been sug- 
gested as an abundance estimator for low metallieity galax- 
ies by Storclii-Bergman, Calzetti & Kinney (1994; see also 
van Zee et al. 1998; Raimann et al. 2000; Denicolo, Ter- 
levicli & Terlevich 2002). However, for abundances greater 
than 0.5 solar, a "fold" in the behavior of [N II]/Ha is 
clearly present in the diagnostic plot of [O IIIJ/H/3 vs. 
[N II]/Ha, as revealed by Dopita et al. (2000). 



puted using the BC03 code. From there, the SED 
corresponding exclusively to MS ionizing stars is 
used to derive r] and Qh° , the number of ionizing 
photons, which then allows us to infer the nebu- 
lar emission-line spectrum by simple interpolation 
within the grid of H II region models, as described 
in §2.4.1. Finally, each emission-line equivalent 
width {t) of line i at wavelength Aj correspond- 
ing to the above evolutionary age t is derived by 
dividing the line luminosity L\{t) by the SED lu- 
minosity L*{Xi, t) of the integrated stellar popula- 
tion (i.e. all stars included) averaged over a fixed 
narrow band-pass A^^ . 

2.4.4- The concept of "equivalent H 11 region" 

In the present model, we are assuming that the 
integrated emission line spectrum of galaxies is 
well represented by an "equivalent H II region," 
ionized by a cluster of MS stars with the same 
mass spectrum as the whole galaxy. In support of 
this assumption, it has been shown by Kobulnicky 
et al. (1999) that the integrated [O II] and [O III] 
line luminosity mimics that of individual H II re- 
gions, and that the chemical abundances of the 
gas phase derived from individual nebulae or from 
integrated galaxy spectra coincide at least within 
the intrinsic uncertainties of the proposed abun- 
dance determination method. For ages older than 
10 Myr and for a continuous star formation rate, 
the number of newly born massive stars relative 
to those having aged and producing less ionizing 
photons, achieve a near steady state regime of con- 
stant hardness, as can be seen in the right-hand 
panel of Figure 1, which shows the time evolution 
of 7] for different continuous SFR. In this adopted 
scenario, the shape of the ionizing spectrum de- 
pends therefore only on the IMF and the metal- 
lieity of the stellar population. We are implicitly 
assuming that the error introduced by calculating 
the nebular spectra using a unique (but variable) 
galaxy age > 10 Myr (rather than taking into ac- 
count a distribution of H II region ages) , is small 
in comparison with the uncertainties introduced 
by the other nebular parameters. Actually, our 
models are intended for comparison with normal 
galaxies, which are massive enough to minimize 
the dispersion because of statistical fluctuations 
in the sampled population. However, as has been 
pointed out by Cervino et al. (2002), this question 
is susceptible to be investigated in more detail. 
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Interestingly, if one seeks the equivalent m„p 
of a zero-age instantaneous burst (Fig. la), which 
would reproduce the same SED hardness rj as that 
from the corresponding asymptotic value in the 
continuous SFR case with m„p = 125M0 (Fig lb) 
and with the same Z, one finds that this zero-age 
equivalent m„p is ~ 60Mq. This result was con- 
firmed for all the metallicities studied in this pa- 
per. It follows that the continuous SFR models of 
Fig. 3b with rriup = 125M0 are, to all practical 
purposes, equivalent to a zero-age instantaneous 
burst with m„p ~ 6OM0. 

3. Model results and discussion 

We first discuss the principal results of our 
model predictions based on exponentially declin- 
ing SFRs and then proceed to compare these with 
the spectrophotometric atlas of JFFC. 

3.1. An atlas of calculated equivalent 
widths for non-instantaneous bursts 

Table 1 summarizes the contents of our model 
grid ADEMis^. The models listed correspond to 
the Salpeter IMF, with rriup = 125M©, and a SFR 
^{t) exponentially declining with time: 

*(0 = r-iexp(-i/T) , (1) 

with T = 5, and 9 Gyr, as indicated. The gas and 
metallicity of the stellar population are the same 
in all models of Table 1, covering the two cases 
of = Zgas = Z = 1 as well as 0.4 Zq. The 
ionization parameter is log U = —3.00. In the 
electronic version, the models cover a wider range 
of the following parameters: r = 1, 3, 5, 7, and 9 
Gyr, and cxd (i. e. constant SFR), log U = —3.25, 
-3.00, -2.75, -2.50 and -2.25, rriup = 40, 60, 
80 and 125 Mq and Z = 0.2, 0.4, 1.0 and 1.5 Zq 
(§2.2). 

In Table 1, column (1) gives the metallicity Z, 
column (2) r (see eq. 1), column (3) m„p (in units 
of Mq), column (4) logC/, column (5) the age t 
of the burst (in units of years), column (6) the 
hardness of the ionizing MS spectrum 77. In col- 
umn (7) , we list Qh° 7 the ionizing photon lumi- 
nosity of MS stars (s~^) per unit solar mass of 
the parent galaxy total mass, column (8) contains 

^Also available as an ASCII file via ftp to 
http: / / www.cida.ve /~ magris/ademis . 



the absolute visual magnitude, My, per unit solar 
mass of the parent galaxy mass. In columns (9)- 
(13) we present the equivalent widths of [O II], H/3, 
[O III], Ha, and [N II] (sec §2.4.3). The emission- 
line fluxes of [Oil], H/3, [O III], Ha, and [Nil], 
shown in columns (14)-(18), are in ergs s~^ per 
unit solar mass of the parent galaxy mass. 

3.2. Relationships between EWs cind metal- 
licity Z or age t 

The relationship between W''a([0 II]) and I'V'A(Ha) 
as a function of time for an exponentially declin- 
ing SFR (Eq. 1) with a fixed r = 3 Gyr is shown 
in Fig 4. Each line joins models of increasing 
age t, from right to left, keeping the ionization 
parameter constant. The three curves in each in- 
set correspond to different ionization parameters, 
covering the values: log U = —2.75, —3.00 and 
—3.25 (long-dashed curve). Models for four dif- 
ferent upper mass limits m„p of 125, 80, 60 and 
40 Mq, are depicted in the four panels of Figs. 
4a-d, respectively. The model metallicity varies 
from inset to inset, in a vertical fashion, to cover 
the four cases Z = 0.2 (highest inset), 0.4, 1.0, 
and 1.5 (lowest inset). As a reference, a star and 
a circle are used to indicate ages of 1 and 12 Gyr, 
respectively''. Filled squares represent the JFFC 
sample (§2.3). The same models are shown in Fig- 
ures 5a-5d {rriup) for the case of Wa([0 III]) versus 
Wxi^a), while for VFa([N II]) versus VFa([0 II]), 
they appear in Figs. 6a-6d. In Figure 7, we plot 
Wa([0 III])-hWA([0 II]) versus WA(Ha), assuming 
the value to„j, = 125Mq. We adopt a logarith- 
mic scale in our figures since Jansen et al. (2000b) 
pointed out that their EWs plots show a near 
uniform dispersion when expressed in logarithmic 
units. 

It is interesting to analyze Figs 4-7 in the light 
of the discussion of the data presented in §4.2 of 
Jansen etal. (2000b). We concluded the follow- 
ing^: 

'^The model ages can be read from Figure 9 (curves labeled 
T = 3), which displays the Ha EWs as a function of t (see 
§3.4). 

®In many line ratio diagnostic diagrams in the literature, 
the H/3 or Ha line appears in the denominator. In our EW 
plots involving Wx(Ha) (Figs. 4, 5, 7 and 8), a normaliza- 
tion is also occurring implicitly. In the case of line ratios, 

the luminosity of the stellar burst (via Qjjo) is taken out 
when dividing by H/3. To some extent, this occurs when 
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1. Because high metaUieity models appear too 
far below the data in Figs 4a-d when m„j, < 
8OM0, we consider that our best case contin- 
uous SFR model favors high values for m„p. 
For the sake of the discussion, we will as- 
sume m„p = I25M0 for the remainder of 
this Section and the next. 

2. Although we cannot infer in a simple way 
the metallicity using EW diagrams, a com- 
parison of models of different metallicities in 
Figure 4a and 5a tells us that galaxies with 
a high Wx{[0 II]) and/or a high Wx{[0 III]) 
(i.e. objects toward the upper right corner) 
could not be fitted satisfactorily using solar 
or above metallicity models. Therefore, our 
models suggest abundances lower than solar 
for this subset. This indication is consistent 
with the findings of Jansen et al. (2000b), 
who pointed out that these high iyA([0 II]) 
objects in fact correspond to a lower lumi- 
nosity (and lower metallicity) subset of their 
survey^. 

3. The data show a trend of increasing W\ ( [O II] ) 
or Wa([OIII]) with Wx{Ra). Our models 
suggest that these could in part be due to 
a younger age effect, since the solid lines 
present an angle of inclination that is com- 
parable to the trend in the data. However, 
the increase in Wx{[0 II]) or Wx{[0 III]) [but 
not that of Wa(Hq:)] may be caused as well 
by a systematic trend in metallicity, since 
models shift vertically by large amounts 
when varying Z, as shown by the insets of 
Fig 4a. This is not surprising, given that the 
luminosities of coUisionally excited lines are 

Wx(iioi) (or any emission-line equivalent width) is plotted, 
since it implies a division of the line luminosity by the un- 
derlying continuum luminosity. A new variable is implicitly 
introduced in the models, however, namely the burst age t 
(or equivalently r), because there exist two major contri- 
butions to the underlying stellar continuum: the youngest 
hot stars, which scale with Qjjo , and the older stellar pop- 
ulation, which does not synchronize with the temporal evo- 
lution of the emission lines. 

Jansen ctal. (2000b) established convincingly that the 
galaxies with a [O II] vs. Hq+[N II] relation steeper (in lin- 
ear plots) than the mean value of 0.4 previously reported by 
Kcnnicutt (1992) corresponded in fact to a lower metallic- 
ity galaxy subset if one adopts the calibration of Zaritsky, 
Kcnnicutt, & Huchra (1994), and Kobulnicky & Zaritsky 
(1999) to infer the galaxy metal abundances. 



strongly dependent on the gas temperature, 
which, in turn, is mostly a function of gas 
metallicity. We conclude that variations in 
age t and in metallicity Z are the two most 
significant contributors to the trends seen in 
the data. 

4. Even though a variation in the age t of 
models mimics part of the trend observed 
in the data in Fig 4a, we cannot reli- 
ably estimate the age (or the metallicity) 
using Wa([0 II]). For instance, galaxies 
in the Jansen et al. (2000b) sample with 
Wx{[0 II]) > 30 A are in general expected to 
possess subsolar abundances with Z < 0.5 
Zq. [In that case, their observed WA(Ha) 
appears confined to values below ^ 70 
A.] However, much larger EW([0 II]) than 
present in the current sample are in prin- 
ciple possible, but only if the Ha equiva- 
lent widths were correspondingly very large. 
In our scheme, such an object would cor- 
respond to stellar bursts characterized by 
a very young stellar population. The pre- 
dicted Wx(H.a) in that case can be found 
amongst the models of Table 1. 

5. In the case of the Ha recombination line, 
as expected, its EW is independent of U, 
while (by definition) high- (low-) U models 
are characterized by higher (lower) [O III] 
EWs and lower (higher) [O II] EWs. 

6. In models, WA([Oni]) is much more sen- 
sitive to variations in U (compare Fig. 5a 
with 4a) than Wa([0 II]). This is consistent 
with the range of values taken by Wa ([O III]) 
within the data, which are much larger than 
that of H^a([0 II]). We recah that the hard- 
ness T] becomes constant after an age as short 
as w 10^ yr, as shown in §2.4.1, and there- 
fore, changes in Wa([0 III]) are not due to 
changes in r], but are entirely due to changes 
in Qh'i and in the integrated stellar contin- 
uum level. 

7. The [N II]-[0 II] EW plot of Fig. 6 indicate 
that variations in age t run almost vertically 
and are therefore nearly orthogonal with re- 
spect to variations in metallicity Z. Further- 
more, the models are not sensitive to the ion- 
ization parameter. This EW diagram might 
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therefore serve as an abundance diagnostic, 
a property worth investigating further. 

8. The dependence on the ionization parameter 
is almost absent in the plot of VF;v([0 III]) + 
Wa([0 II]) against Wx{lia), as shown in Fig. 
7. 

3.3. The abundance— m„p degeneracy 

In order to exemplify the degeneracy between 
abundance and m„p (see Sec. 2.4.1) in the case 
of continuous star formation, we present in Figure 
8 the loci occupied by models with metallicities 
0.4, 1.0, and 1.5 times solar, in a plane defined by 
(panel a) VFa([0 II]) versus VFA(Ha) and (panel 
b) Wa([OIII]) versus WxiHa). Areas with dif- 
ferent shades correspond to models with different 

values (as labeled). The other parameters are 
log U between —3.25, and —2.50, a SFR declining 
exponentially with time (Eq. 1) with r anywhere 
between 1 Gyr, and oo (i.e. constant SFR), and 
a Salpcter IMF with m„p= 125, 80, and 60 Mq. 
These parameters encompass a wide range of pos- 
sibilities for the galaxies considered. The main is- 
sue to highlight here is the extent to which metal- 
licity and m„p are degenerate as showed by the 
significant overlap between the different shaded ar- 
eas in the [O III] and Ha EW plot of Figure 8b. 
The abundance-m„)5 degeneracy^" must therefore 
be taken into account for the analysis of emission 
line properties, especially in small [O III] equiv- 
alent widths galaxies, in order to determine the 
origin of the low excitation. Because of the low 
ionization potential of O to 0+, [O II] does not 
suffer from such a degeneracy and it is therefore 
a better indicator of abundances in emission-line 
galaxies. 

The comparison of our models with observed 
line ratios allows us to conclude that the metallic- 
ities of the galaxies in the JFFC sample do not ex- 
tend beyond 1.5 ^Qand that an important fraction 
of them actually belongs to a galaxy population of 
lower metallicity than solar. 

'The degeneracy in the case of continuous SFR is at least 
partly multivariate since the new vajriable t (via Eq. 1) 
malies the SED (and hence the EWs) a function of time 
due to the evolution of the stellar continuum (even though 
the hardness rj is constant after lO'^ yr) . We did not explore 
in detail the boundaries of the t-Z^ruup degeneracy. 



3.4. Varying the rate of SFR decline r 

Models with different e-folding time r (but 
equal Z and U) show relationships between the 
[O II] or the [O III] EWs, and the Ha EWs, very 
similar to those studied in §3.2, except for their 
temporal evolution rate, which scales differently. 
This is expected since the hardness of the ioniz- 
ing SED is constant after w 10^ yr, as shown in 
§2.4.1, and models of interest are typically much 
older than this. Also, since our EW models corre- 
spond to the properties of the whole host galaxy, 
the SED in the optical, against which the EW are 
defined, is not sensitive to the youngest star colors. 

The equivalent width-age relationship, as ex- 
pected and mentioned above, strongly depends on 
the adopted particular SFR history, as illustrated 
by Fig. 9, in which we plot the time evolution of 
Wx (Ha) for similar models than shown in previous 
figures except that r is now varied. We only con- 
sider the case niup = 125 M©, which is the model 
favored by our discussion in §3.2. The horizon- 
tal lines circumscribe the region occupied by the 
JFFC galaxies. The WA(Ha) versus age relation 
is not sensitive to the ionization parameter. How- 
ever, it shows a residual dependence on metallicity, 
as can be appreciated in Figure 9. This is simply 
due to variations in the mean effective temper- 
ature of the ionizing stars with Z^, which brings 
about small changes in the number of ionizing pho- 
tons and also to variations in the integrated stellar 
continuum. 

The first conclusion reached from introspection 
of Fig. 9 is that even for the most active galax- 
ies (in the sense of higher star- forming rates), the 
stellar population must be older than w 1 Gyr in 
order to fit within the observed EW range. Actu- 
ally, even allowing for a very long-lasting episode 
of high SFR, i.e. models with high T-valucs, the 
requirement that these models lie below the ob- 
served 100 A T40\(Ha) upper limit, imply such a 
lower limit on the burst age, as expected in a sam- 
ple of normal field galaxies. On the other hand, 
we would conclude that galaxies from other sam- 
ples, which would be observed to have WA(Ha) 
> 100 A. must therefore be actively forming stars, 
and their dominant stellar population could not 
be much older than 1 Gyr. 

The second conclusion derived from Fig. 9 is 
that galaxies found with a small VFA(Ha) (< lOA) 
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arc currently forming stars at a very slow rate with 
respect to their previous SFR. This is implied by 
the requirement of a small r (< 3 Gyr) so that the 
models can reproduce M^A(Ha) < 10 A. 

A more precise determination of the star forma- 
tion history is not possible, using only Vt^A(Ha). 
As is obvious in Fig. 9, the same WA(Ha) is ob- 
tained by either increasing or decreasing t and r 
simultaneously. The degeneration between both 
parameters can be overcome by using specific trac- 
ers. Kauffmann et al (2003), for instance, used the 
strength of the 4000 A break as well as the stellar 
absorption line index H(5 to constrain the star for- 
mation history of the galaxies. Recently, J. Mateu, 
G. Bruzual A., & G. Magris C. (2003, in prepa- 
ration) have model-fitted the stellar continuum of 
the JFFC sample and were able to infer the tempo- 
ral evolution of the SFR of the underlying stellar 
population for most of the sample. 

As for the line ratios [O II]/H/3, [O III]/H/? and 
([O 11] + [0 III]/)H/3, their behavior as a function 
of the Ha EW is shown in Figure 10 for models 
with a Salpeter IMF and niup = 125 Mq. Each 
strip represents the loci of models of the same 
metallicity with log U anywhere between —3.25 
and —2.25 and t between 1.0 and infinity (equiv- 
alent to a constant SFR). The corresponding ages 
of the models can be inferred from Fig. 9. Since 
models with Z — 0.2 and 0.4 would overlap, we 
only show the case Z = 0.4. We infer from panel a 
that the [O Il]/H/j ratio could give us a reasonable 
estimate of abundances, specially for metallicities 
above solar, while the classical abundance indica- 
tor ([O II]-F[0 III])/H/3 ratio (Pagel et al. 1979), 
as expected, remains a better diagnostic since the 
shaded strips are much thinner in the correspond- 
ing panel c. The [O III]/H/3 ratio in panel b, re- 
mains a poor estimator of metal abundance in the 
nebulae. 

3.5. General discussion 

Stasinska et al. (2001) presented instantaneous 
burst models for the evolution of H II galaxies 
and suggested that the presence of an underly- 
ing old stellar component was required to lower 
the Wa(H/3) predicted by their models sufficiently. 
Similarly, Moy et al. (2001) have shown that an 
instantaneous or a very short-lasting burst of star 
formation predicts very high hydrogen lines equiv- 
alent widths. Therefore, they added an underly- 



ing population to their starburst galaxy model. 
Our integrated continuous SFR model includes 
and evolves its own intrinsic older stellar popu- 
lation in a self-consistent manner. This has the 
effect of "raising" the continuum in the H/3 and 
Ha wavelength regions, thereby lowering the neb- 
ular emission line equivalent widths. 

An important issue to emphasize is that the 
ADEMis model is able to simultaneously account 
for the observed oxygcn-to-hydrogen emission-line 
ratios as well as for the Ha equivalent widths. In 
any event, the presence of an older population af- 
fects EWs but not the line ratios"'^"'^. As has been 
explained in §2.4.1, in the case of a continuous 
SFR, the hardness of the ionizing spectrum evolves 
only during the first few Myr, until a balance be- 
tween stellar births and deaths is established. In 
this latter regime, the parameters that govern the 
emission line ratios are the metallicity and the 
IMF upper mass limit, and, to a lesser extent, the 
slope of the initial mass function (a parameter not 
explored in this paper). 

Even though our approach is not identical to 
that of Chariot & Longhetti (2001), both models 
have in common that they decouple the evolution 
of the ionizing stellar population from the evolu- 
tion of the older stellar population since they re- 
quire us to select, for a fixed metallicity, the hard- 
ness of ionizing spectrum. This is done in our 
case by appropriately selecting the most suitable 
IMF upper mass limit, or, in the case of Chariot 
& Longhetti, by selecting the stellar burst age. 

4. Conclusions 

We have presented a new library of models, 
ADEMIS, which combine the stellar population his- 
tory with the emission line-strengths. These mod- 
els can be used to interpret the properties of ob- 
jects that arc characterized by either a very recent 
and strong burst or by a secular and lower star for- 
mation rate. 

Starting from simple assumptions about the 
distribution of H II regions in disk galaxies, 
namely, a uniform distribution of dust-free gas 
and stars, in which young stars arc distributed in 
mass in a similar fashion as the older stars, we 

The line ratios as a function of age depend only on rj, since 
both U and Z are assumed constant with time. 
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have shown how it is possible to compute several 
properties of star forming galaxies. A photoion- 
ization code is used to compute the intensity of the 
emission lines, assuming no temporal evolution of 
U, while an evolutionary synthesis code provides 
the number luminosity and hardness of the MS 
ionizing photons. These, in conjunction with the 
photoionization code and the integrated stellar 
distribution, allow us to compute the temporal 
evolution of the emission line equivalent widths. 

Diagnostic diagrams built by comparing pairs 
of equivalent widths of various lines, show that our 
model predictions fall in figures in the same region 
as the observed galaxies of JFFC. In particular, we 
conclude the following. 

1. A factor of 10 in [/ and a range in metallicity 
between 0.2 and 1.5 Zq seems appropriate to 

describe the observed properties of a repre- 
sentative sample of nearby galaxies. 

2. The hardness of the ionizing continuum is 
constant after 10^ yr in case of a continuous 
SFR. 

3. Our favored continuous SFR model with 
m„p = 125M0 presents a hardness r] very 
similar to that produced by an instanta- 
neous burst with an IMF having only m„p = 

6OM0. 

4. Galaxies with VFa([0 II]) > 30 A and 
lyA(Ha) ^ 70 A appear to have subsolar 
abundances. 

5. Galaxies with WxCHa) < 100 A have a dom- 
inant underlying stellar population older 
than at least 1 Gyr. Only galaxies with 
a recent burst of star formation {t < 1 Gyr) 
are able to show Ha EW in excess of 100 A. 

6. The information provided by the [O III] line 
is of limited use, owing to the strong niup-Z 
degeneracy, which complicates the interpre- 
tation of this high-excitation indicator, par- 
ticularly in the case of low [O III] EWs. 
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Fig. 1. — (a) Behavior of ry, the hardness of the 
ionizing continuum (as defined in §2.4.1), as a 
function of niup, assuming an instantaneous burst 
with a Salpeter IMF and one of the following 
metallicities: = 0.2, 0.4 , 1.0 and 1.5 Zq. (h) 
Temporal evolution of the hardness for the same 
range in metallicity (as labeled) in the case of 
a continuous SFR decreasing exponentially with 
time, a Salpeter IMF and the same m.up of 125 
Mq. After 10 Myr, i) reaches an asymptotic value, 
which differs depending on metallicity, but can be 
shown to be independent, of the SFR provided 
that the SFR is continuous and either constant or 
monotonically decreasing. The dashed line corre- 
sponds to an instantaneous burst of solar metal- 
licity with the same IMF and m„p. 



Fig. 2.— Diagnostic plot of log ([O III]/H/?) 
against log ([O II]/H/3) for (a) Z = 0.2 Zq, (h) 
Z = 0.4 Zq, (c) Z = 1.0 Zq, and (d) Z = 1.5 Zq. 
The theoretical grid models were calculated with 
an ionizing spectrum from an instantaneous burst 
of age zero, a Salpeter IMF, and niup = 125, 100, 
90, 80, 70 and 60 M©, respectively. The models 
move up in the diagram with increasing m„p [as 
labeled in panel (b)]. Solid lines join models {en- 
circled dots) of equal m„j,, i.e. constant hardness, 
while dotted lines join models of equal ionization 
parameter U. From left to right, log U takes on 
the values of -2.25, -2.50, -2.75, -3.00, -3.25, 
respectively. Filled squares represent the observed 
reddening corrected ratios from the JFFC atlas 
(see §2.3). The arrow in (a) represents the dered- 
dening vector corresponding to E(B-V) = 0.2. 
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Fig. 3.— Diagnostic plot of log([0 III]/[0 II]) 
against log ([N II]/[0 II]) (a): for an instantaneous 
burst at zero age, and (h): for a continuous SFR 
at age 10 Myr, at which point the ionizing hard- 
ness T] has reached its asymptotic value (see Fig. 
lb). In both panels, a Salpeter IMF with the same 
rriup = 125M0 was assumed, and the nitrogen 
abundance was increased via secondary nucleosyn- 
thesis, as explained in §2.2. In the grid shown in 
each panel, models are joined by solid lines with 
the metallicity increasing from left to right, while 
the ionization parameter (log U) increases going 
up. Filled squares represent the observed redden- 
ing corrected ratios from the JFFC atlas. The ar- 
row in (a) represents the dereddening vector cor- 
responding to E(B-V) = 0.2. 



Fig. 4. — Relation between the equivalent widths 
of T4^a([0 II]) and VFa(Hq;). Lines join models of 
increasing age t, from right to left, that were cal- 
culated assuming a continuous exponentially de- 
creasing SFR (Eq. 1) assuming t = 3 Gyr. The 
upper mass limit of the Salpeter IMF varies from 
panel to panel as follows: (a) ni^p = 125 M0, 
(b) ruup = 80 M0, (c) rriup = 60 Mq and (d) 
mup = 40 Mq. Within each panel, different in- 
sets correspond to different mctallicities Z of mod- 
els, increasing from top to bottom: Z — 0.2 Zq, 
Z = 0.4 Zq, Z = 1.0 Zq, and Z = 1.5 Zq. As a 
guide, an asterisk and a circle indicate the position 
of models with ages of 1 and 12 Gyr, respectively. 
The three lines correspond to a different ioniza- 
tion parameter each and cover the values log U = 
—2.75, —3.00 and —3.25. The long-dashed curve 
represents the lowest U-value model. The filled 
squares represent reddening corrected ratios from 
the JFFC atlas, who reported ±3% to 10% for 
strong emission lines, and ±30% for W\ smaller 
than 2 A. 
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W,(Ha) (A) 

Fig. 6.— Same as Fig. 4, but for VKa([NII] A 
Fig. 5. Same as Fig. 4, but for M^a([0 III] A 6584) vs. Wa([0 II]). The nitrogen abundance 

5007) vs. Wx(Ha). The long-dashed line repre- was increased via secondary nucleosynthesis, as 

sents the lowest U-value model. explained in §2.2. The long-dashed line represents 

the lowest U-value model. 



16 




0.1 

1 10 100 

W^(Ha) (A) 



Fig. 7.— Same as Fig. 4, but for 
T4^a([0 III])+VFa([0 II]) vs. Wx{[0 II] AA 3727) as- 
suming a Salpeter IMF with m„p = 125M0. 
The long-dashed hne represents the lowest U- value 
model. 




w^(Ha) (A) 

Fig. 8.— Relation between (a) T/Fa([0 II] AA 
3727) vs. WxiHa) and (b) Wa([0 III] A 5007) vs. 
Wx(tla). Each shade corresponds to the loci of 
continuous SFR models with the same metallicity 
(as labeled), but with log U comprised between 
—3.25 and —2.50, r ranging between 1 Gyr, and 
oo, and m„p taking the values of 125, 80, and 60 
M.Q. The age t increases from right to left. The 
three metallicities covered are Z = 0.4, 1.0 and 1.5 
Zq . The values observed by JFFC are represented 
in each panel by filled squares. 
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w,(H<«) (A) 



Fig. 9. — Evolution of Wx (Ha) with age t for mod- 
els with a SFR declining exponentially with time, 
a Salpeter IMF, and niup = 125 Mq. The rate 
of SFR decline is represented by r, which takes 
on the values r = 1, 3, 5, and 9 Gyr and oo (i.e. 
constant SFR), as labeled. The dotted line cor- 
responds to models with Z = 0.2 Zq, the short 
dashed line to Z = 0.4 Zq, the solid line to Z = 1.0 
Zq, and the long dashed line to Z = 1.5 Zq. The 
shaded area represents range of 97% of observed 
(HQ!)-values from the JFFC sample of galax- 
ies [three objects exceed the 100 A bar, having a 
Wxlua) of 169, 288 and 328 A]. 



Fig. 10. — Relation between the line ratios (a) 
[O II]/H/3, (b) [O III]/H/?, (c) ([O II] + [0 III])/H/3, 
and the Ha equivalent width. In each panel, 
strips with different shades correspond to models 
with a different metallicity Z, as labeled. Each 
shade indicates the loci of models with an expo- 
nentially declining SFR between r = 1 Gyr and 
00. All models assume a Salpeter IMF and the 
same m„p = 125M0. The relative thickness of 
each strip is mostly due to the dependence on the 
ionization parameter, which here spans the range 
—3.25 < log U < —2.50. The age t increases from 
right to left and can be read from Fig. 9 for any 
particular value of r. Filled squares represent the 
JFFC sample. 
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Table 1 

Model Emission-Line Equivalent Widths and Fluxes 



Emission — Line Equivalent Widtlis (A) log Emission — Line Fluxes (erg s ^^^0 ^) 

Age T, log QjjQ My^ [O II] H(3 [O III] Ho [N II] [O II] H/3 [O III] Ho [N II] 

(Gyr) s~'^Mq~^ A3727 A6007 A6583 X3727 A5007 :X65S3 



x = 5 Gyr Z = 0.4Zq log U = -3.00 



0.5 


X 10-3 





.31 


42 


.49 


11.45 


1713. 


96 


812. 


99 


3030. 


13 


7049, 


33 


1229, 


11 


30 


,81 


30 


.13 


30 


.65 


30 


.58 


29.83 


1.0 


X 10-3 





.31 


42 


.79 


10.70 


1707. 


53 


810. 


10 


2987. 


,53 


7024. 


,48 


1224. 


.76 


31 


.11 


30. 


.42 


30 


.94 


30 


.88 


30.13 


0.5 


X 10-^ 





.22 


43 


.32 


8.46 


783. 


.72 


352. 


06 


757. 


18 


2960. 


,98 


516. 


.38 


31 


.63 


30. 


.95 


31 


.24 


31 


.41 


30.66 


0.1 


X 10-^ 





.21 


43 


.34 


7.54 


468 


.70 


171. 


76 


332. 


61 


1154. 


,30 


201. 


.92 


31 


.66 


30. 


.98 


31 


.22 


31 


.43 


30.68 


0.5 


xio-^ 





.21 


43 


.34 


6.67 


250. 


.54 


78 


.35 


149. 


24 


486. 


,97 


85. 


,22 


31 


.65 


30 


.98 


31 


.22 


31 


.43 


30.68 


0.1 


xio" 





.21 


43, 


.33 


6.36 


206 


.15 


58, 


.68 


111. 


53 


358 


92 


62. 


81 


31 


.65 


30 


.97 


31 


.21 


31 


.43 


30.67 


0.5 


x lo" 





.21 


43 


30 


.'5.70 


139 


53 


30 


67 


57. 


67 


173 


08 


30 


29 


31 


.61 


30 


,94 


31 


,18 


31 


.39 


30.64 


1.0 


X lo" 





.21 


43 


25 


5.51 


121 


50 


23 


.84 


44. 


53 


126 


85 


22 


20 


31 


.57 


30 


,89 


31 


,13 


31 


.35 


30.59 


2.0 


X 10° 





.21 


43, 


.17 


5.41 


105 


.32 


18, 


.54 


34. 


40 


88. 


27 


15. 


,45 


31 


.49 


30 


,81 


31 


,05 


31 


.26 


30.51 


4.0 


X10° 





.21 


42 


.99 


5.51 


90 


.32 


14 


.28 


26. 


34 


60. 


50 


10. 


,59 


31 


.31 


30 


,63 


30 


.88 


31 


.09 


30.33 


6.0 


xio" 





.21 


42 


.82 


5.70 


80 


,11 


11 


.70 


21. 


50 


45. 


93 


8. 


04 


31 


.14 


30 


,46 


30 


.70 


30 


.92 


30.16 


8.0 


X lo" 





.21 


42 


65 


5.91 


71 


,70 


9 


79 


17, 


95 


36. 


28 


6 


35 


30 


.96 


30 


,28 


30 


.53 


30 


.74 


29.99 


1.0 


X lol 





,21 


42 


47 


6.13 


63 


81 


8, 


16 


14 


93 


28. 


74 


5 


03 


30 


.79 


30 


,11 


30 


.35 


30 


.57 


29.81 


1.2 


X 10^ 





.21 


42 


.30 


6.33 


56 


17 


6 


73 


12 


29 


22, 


61 


3. 


96 


30 


.62 


29 


,94 


30 


.18 


30 


.40 


29.64 


1.4 


XIO^ 





.21 


42 


.12 


6.52 


48, 


79 




47 


9 


98 


17 


63 


3. 


,09 


30 


.44 


29 


,76 


30 


.01 


30 


.22 


29.47 



T = 5 Gyr Z =Zq log U = -3.00 



0.5 


X 10^-^ 


0, 


22 


42 


51 


11.29 


632 


49 


722. 


13 


484, 


56 


6671 


10 


2513 


29 


30 


,44 


30 


14 


29. 


92 


30 


,62 


30.19 


1.0 


X 10-3 


0. 


22 


42 


,80 


10.51 


596, 


,35 


696. 


,85 


436. 


64 


6430, 


04 


2398. 


27 


30 


,73 


30 


,44 


30. 


19 


30 


.91 


30.48 


0.5 


X 10-^ 


0. 


17 


43 


,25 


8.21 


200, 


,84 


243. 


36 


82. 


78 


2139. 


.10 


741. 


.40 


31 


,12 


30 


.89 


30. 


38 


31 


.37 


30.91 


0.1 


X 10-^ 


0. 


16 


43, 


,26 


7.49 


123 


,64 


136. 


00 


42. 


22 


988 


10 


339. 


00 


31 


,12 


30 


.91 


30. 


35 


31 


.38 


30.92 


0.5 


X 10-1 


0, 


16 


43 


26 


6.72 


74 


85 


69 


52 


21. 


21 


437 


01 


149 


74 


31 


,12 


30 


,91 


30. 


35 


31 


.38 


30.91 


0.1 


X lo" 


0, 


16 


43 


26 


6.42 


62 


55 


52 


77 


16, 


05 


324 


48 


111 


18 


31 


,12 


30 


,90 


30, 


34 


31 


.37 


30.91 


0.5 


X lo" 


0. 


16 


43 


,22 


5.80 


44 


,37 


28 


,30 


8 


,53 


164. 


38 


56. 


32 


31 


,08 


30 


,86 


30, 


,31 


31 


.34 


30.87 


1.0 


xio" 


0. 


16 


43 


,18 


5.63 


39 


,32 


22 


,22 


6 


,66 


121. 


26 


41. 


.55 


31 


,04 


30 


,82 


30. 


,26 


31 


.30 


30.83 


2.0 


xio" 


0. 


16 


43 


,09 


5.55 


34, 


,85 


17, 


,70 


5 


,28 


87. 


03 


29. 


82 


30 


,95 


30 


,74 


30. 


18 


31 


.21 


30.74 


4.0 


X lo" 


0, 


16 


42 


92 


5.69 


30 


83 


14 


04 


4 


17 


61. 


56 


21 


09 


30 


,78 


30 


,56 


30 


00 


31 


.03 


30.57 


6.0 


X lo" 





16 


42 


75 


5.90 


28 


06 


11 


76 


3 


48 


47, 


67 


16 


33 


30 


,60 


30 


,39 


29 


83 


30 


.86 


30.40 


8.0 


X lO" 


0, 


16 


42 


,57 


6.13 


25 


,66 


10, 


,01 


2 


,96 


38, 


23 


13 


10 


30 


,43 


30 


,21 


29 


66 


30 


.69 


30.22 


1.0 


XIO^ 


0. 


16 


42 


.40 


6.36 


23 


.35 


8 


,47 


2 


,50 


30. 


64 


10. 


.50 


30 


,26 


30 


,04 


29. 


,48 


30 


.51 


30.05 


1.2 


XIO^ 


0. 


16 


42 


.22 


6.58 


21 


.08 


7, 


,09 


2 


,09 


24. 


38 


8. 


.35 


30 


,08 


29 


.87 


29. 


,31 


30 


.34 


29.88 


1.4 


XIO^ 


0. 


16 


42 


.05 


6.78 


18 


.81 


5, 


,85 


1 


,73 


19. 


16 


6. 


.57 


29 


,91 


29 


.69 


29. 


,14 


30 


.17 


29.70 



















T = 9 Gyi 




Z = 


0.4ZQ 


log 17 = - 


3.00 




















0.5 


X 10- 


3 





.31 


42 


.23 


12.09 


1713.96 


812. 


99 


3030.13 


7049.33 


1229. 


,11 


30 


.55 


29 


.87 


30.40 


30 


.33 


29.57 


1.0 


X 10- 


3 





.31 


42 


.53 


11.33 


1707.53 


810. 


10 


2987.53 


7024.48 


1224. 


76 


30 


.85 


30 


.17 


30.69 


30 


.63 


29.87 


0.5 


X 10- 


2 





.22 


43 


06 


9.09 


783.78 


352. 


10 


757.31 


2961.28 


516 


43 


31 


,38 


30 


70 


30,99 


31 


.16 


30.40 


0.1 


X 10^ 


1 





.21 


43 


08 


8.18 


468.82 


171. 


81 


332.74 


1154.68 


201 


99 


31 


,40 


30 


72 


30,97 


31 


.18 


30.42 


0.5 


X 10- 


1 





.21 


43, 


.08 


7.30 


250.82 


78, 


,45 


149.46 


487.69 


85. 


35 


31 


,40 


30 


.72 


30.97 


31 


.18 


30.42 


0.1 


xio" 







.21 


43 


.08 


6.99 


206.58 


58 


,84 


111.84 


359.93 


62. 


,99 


31 


.40 


30 


.72 


30.96 


31 


.18 


30.42 


0.5 


xio" 







.21 


43 


.06 


6.31 


140.73 


31 


,06 


58.42 


175.48 


30. 


71 


31 


.38 


30 


.70 


30.94 


31 


.16 


30.40 


1.0 


X lo" 







.21 


43 


04 


6.08 


123.42 


24 


,42 


45.64 


130.39 


22 


82 


31 


.35 


30 


.68 


30.92 


31 


.13 


30.38 


2.0 


X lO" 







.21 


42 


99 


5.91 


108.53 


19 


43 


36.08 


93.52 


16 


37 


31 


.31 


30 


.63 


30.87 


31 


.09 


30.33 


4.0 


X lo" 







.21 


42 


.89 


5.87 


95.89 


15 


,66 


28.94 


68.09 


11. 


92 


31 


.21 


30 


.53 


30.77 


30 


.99 


30.23 


6.0 


xio" 







.21 


42 


.79 


5.93 


88.22 


13 


,56 


24.97 


55.40 


9. 


,69 


31 


.11 


30 


.43 


30.68 


30 


.89 


30.14 


8.0 


xio" 







.21 


42 


.70 


6.03 


82.50 


12, 


,09 


22.23 


47.26 


8. 


,27 


31 


.02 


30 


.34 


30.58 


30 


.80 


30.04 


1.0 


XIO^ 







.21 


42. 


.60 


6.14 


77.55 


10 


,89 


19.99 


40.99 


7. 


,17 


30 


.92 


30 


.24 


30.49 


30 


.70 


29.94 


1.2 


XIO^ 







.21 


42 


.51 


6.26 


73.02 


9 


,84 


18.04 


35.79 


6. 


,26 


30. 


.82 


30 


.15 


30.39 


30 


.60 


29.85 


1.4 


XIO^ 







.21 


42 


.41 


6.38 


68.77 


8 


,90 


16.31 


31.39 


5. 


,49 


30 


.73 


30 


.05 


30.29 


30 


.51 


29.75 



T = 9 Gyr Z =Zq log U = -3.00 



0.5 


X 10-3 


0. 


22 


42 


25 


11.92 


632 


49 


722. 


13 


484. 


56 


6671 


10 


2513 


29 


30 


.19 


29 


89 


29. 


67 


30 


,36 


29.94 


1.0 


X 10^3 


0, 


22 


42 


54 


11.15 


596 


35 


696, 


86 


436 


65 


6430 


05 


2398 


27 


30 


,47 


30 


18 


29, 


93 


30 


,66 


30.23 


0.5 


X 10-2 


0, 


17 


43 


.00 


8.85 


200 


,86 


243, 


39 


82. 


80 


2139 


35 


741, 


49 


30 


,86 


30 


.64 


30. 


12 


31 


,11 


30.65 


0.1 


X IQ-l 


0. 


16 


43 


.01 


8.13 


123 


,67 


136. 


03 


42. 


23 


988. 


,42 


339. 


,12 


30 


.87 


30 


.65 


30. 


10 


31 


.13 


30.66 


0.5 


xio-1 


0. 


16 


43 


.01 


7.36 


74, 


,92 


69 


.60 


21. 


24 


437. 


63 


149. 


95 


30 


.87 


30 


.65 


30. 


09 


31 


.12 


30.66 


0.1 


X lO" 


0. 


16 


43 


01 


7.05 


62 


,67 


52 


91 


16. 


09 


325 


37 


Ill 


48 


30 


,86 


30 


,65 


30, 


09 


31 


.12 


30.66 


0.5 


X lo" 


0, 


16 


42 


99 


6.41 


44 


72 


28 


65 


8 


64 


166 


57 


57 


07 


30 


,84 


30 


,63 


30 


07 


31 


,10 


30.64 


1.0 


X lO" 


0, 


16 


42 


96 


6.19 


39 


88 


22 


75 


6 


82 


124. 


56 


42 


68 


30 


,82 


30 


,60 


30 


05 


31 


,08 


30.61 


2.0 


X lO" 


0, 


16 


42 


92 


6.05 


35 


,79 


18 


.50 


5 


.52 


91. 


97 


31. 


,51 


30 


,77 


30 


,56 


30 


00 


31 


,03 


30.57 


4.0 


X lO" 


0, 


16 


42 


.82 


6.04 


32 


,43 


15 


.29 


4 


.55 


68. 


83 


23. 


,58 


30 


,68 


30 


,46 


29 


91 


30 


,93 


30.47 


6.0 


X10° 


0. 


16 


42 


.72 


6.12 


30 


,40 


13 


.47 


4 


.00 


56. 


93 


19. 


,51 


30 


,58 


30 


,36 


29, 


81 


30 


.84 


30.37 


8.0 


xio" 


0. 


16 


42 


.63 


6.23 


28 


.84 


12 


.16 


3 


.60 


49. 


16 


16. 


,85 


30. 


.48 


30 


.27 


29. 


.71 


30 


.74 


30.28 


1.0 


xio^^ 


0. 


16 


42 


.53 


6.35 


27, 


.47 


11 


.07 


3 


.28 


43. 


02 


14. 


,74 


30 


.39 


30 


.17 


29. 


.61 


30 


.65 


30.18 


1.2 


xio^ 


0. 


16 


42 


.43 


6.48 


26 


.22 


10 


.10 


2 


.99 


37. 


89 


12. 


,98 


30 


.29 


30 


.08 


29. 


.52 


30 


.55 


30.08 


1.4 


XIO^ 


0. 


16 


42 


.34 


6.61 


25, 


.04 


9 


.22 


2 


.73 


33. 


45 


11. 


,46 


30 


.19 


29 


.98 


29. 


,42 


30 


.45 


29.99 



^in mag per galaxy mass in units of solar mass 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample 
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